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1.0 SUMMARY 


This document describes the structure, design, and programming details of the QR 
program. Use of the program is described in volume I of this document. QR can be used 
as a standalone program or as a module of a program system called DYLOFLEX (see 
fig. 1) which was developed for NASA under contract NASl-13918 (ref. 1). 

1.1 PURPOSE AND OBJECTIVES 

The QR program ( L224 ) is primarily a tool for applying classical control systems 
analysis and synthesis techniques. To make use of the program, the physical system to 
be analyzed is initially described by a system of simultaneous ordinary linear 
differential equations with constant coefficients. By manipulating the input matrix of 
equations, sets of eigenvalues that represent the denominator (POLES) and numerator 
(ZEROS) of the system can be calculated. From the transfer functions, studies of time 
response, frequency response, and root locus may be performed. In addition to the 
aforementioned analysis techniques, two unsteady aerodynamics analyses options 
(unsteady frequency response and V-g analysis) are available. 

1.2 PROGRAM DESIGN AND STRUCTURE 

The QR program can be divided into five different analyses: eigenvalue rooting, time 
response, frequency/PSD response, root locus, and unsteady aerodynamics. Each major 
computational task is generally built into one overlay of the program. Figure 2 
illustrates card input, external magnetic file input, punched matrix card output and 
plot vector file output. 




Figure 1. - DYLOFLEX Flow Chert 












2.0 PROGRAM DESCRIPTION 


2.1 DISCUSSION 

A fairly complete program description of QR can be found in sections 4.0, 5.0 and 6.0 in 
volume I of this document. Figure 3 shows a block diagram of the overlay structure in 
QR. The last letter of the overlay block name generally signifies what is actually done 
during execution. The M in QRM stands for main executive program (MAIN); T in QRT 
stands for time response analysis; R in QRR stands for root locus; U in QRU stands for 
unsteady aerodynamics. Sections 2.1.1 through 2.1.17 detail the task of each overlay. 

2.1.1 OVERLAY QRM, 0,0 

This overlay is the QR main program and its purpose is to control the tasks that are 
undertaken. The control process is very straight forward. QR reads a data card from the 
INPUT file and determines if a valid command has been input. Assuming a valid 
command has been input, this overlay will do one of two things: (1) branch to another 
subroutine or position in the main program (QR) where additional input data is read or 
analytic steps are performed, (2) call another overlay to perform time response, 
frequency response, root locus plotting, V-g analysis plotting, or unsteady aerodynamics 
options. 

The principal analytic step undertaken in the QRM, 0,0 overlay is eigenvalue rooting. 
The technique used to reduce the polynomial input matrix and calculate eigenvalues is 
discussed in reference 2. In addition, this overlay: (1) controls looping for phase gain 
locus and calculation of open loop zero and poles {section 4 *2, volume I discusses this 
calculation): (2) controls looping for altitude and reduced frequencies when performing 
V-g analysis (section 4.9 of volume I discusses the theory for V-g analysis); and (3) 
calculates the polynomial coefficients associated with a given set of eigenvalues; 
(section 4.4of volume I discusses this calculation). 

I/O devices are directly or indirectly opened in this overlay. They include: 

(1) INPUT file or TAPE5 

(2) OUTPUT file or TAPE6 

(3) PUNCH file or TAPE93 

(4) Interactive data base file or TAPE7 

(5) Matrix storage file or TAPE90 

(6) *Tn-core” matrix file or TAPE91 

(7) Root locus and flutter eigenvalue storage file or TAPE96 


4 



QRIVU,0 

Initialization 


QRM,7,0 

Program 

completion 


QRT,2,0 
Time response 


QRU,6,0 

Unsteady 

frequency 

response 



Figure 3, 




















(8) SC4020 plot output file or TAPE99 

(9) Plot vector file or TAPE92 

(10) READ MATRIX TAPE data files - TAPEl, TAPE2, TAPE3, TAPE4, TAPES, 
TAPEIO, TAPEll, TAPE12, TAPE13, TAPE14, TAPE15, TAPE16, TAPE17, 
TAPE18, TAPE19, TAPE20 

Figure 4 is a macro flow chart for the QRM,0,0 overlay. 

A map of subroutines called in QRM,0,0 follows (table 1). 

2.1.2 OVERLAY QRM,1,0 

This overlay calculates constants that are to be used in QR, prints an identification 
page on the output file (TAPE6), sets up pooled program buffers, and reads the first QR 
data card from the input file (TAPES). 


Table 2. — Subroutines Called by Overlay (QRM, 1,0} 


QRM10 Calls 


t DYLOFLEX library routine 
* CDC system library routine 


BIGGER 

COST" 

KOMSTRt 

EOF" 

MUSER 

TTY* 

CLOCK" 

WA10 

RWA10 

VUSER 

FETADt 

LOCF" 

STRMOVt 

PRINTV - calls 


USERSB 

RA10 


A10R1 

BUFFIL 
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Figure 4. — Macro Flow Chart of Overlay (QRM, 0,0) 
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Table 1. - Subroutines Called by Overlay (QRM,0,0) 

(QRMAO) 

Program OR 


ACCT 

AUTOFL 

BEGINS 


AUTOFL 

REQFL* 

AUTOFL 

BLKFLL 

CLEAR 

RWA10 

WARN 


REQFLt 

REQFLt 


BLKFLL 

CANCEL / FORMPM 
I RPRINT 


CLEAR 

COL 

COLROW RW12 

RWI5 

CRAMER / RWI5 




SECOND* 

STRMOVt 


RUW (see below) 


RUW (see below) 


RUW 


^ COL 
ENDQR 
ISCANt 
KNVRTt 
KOMPARt 
LOCF* 
SCAMB 
^ STRMOVt 


t DYLOFLEX library routine 
* CDC FORTRAN library routine 


AUTOFL 



Table 1. — {continued} 


CREAD CLEAR 

CSAVE 

ENDQR AUTOFL REQFLt 

QQQQ7 

REGRUP 

RUW (see above) 

RWE12 RUW (see above) 

^ RWI5 RUW (see above) 

CSAVE 

ENDQR AUTOFL REQFLt 

FORMPM 

FRESP AUTOFL REQFLt 

FRGRID / RWE1 2 (see above) 

I RWI5 (see above) 

FTPLOT AUTOFL REQFLt 

GAINDC 

GAINS < RUW (see above) 

RWE12 (see above) 

RWI 5 (see above) 

GUST / RWE12 (see above) 

I RWI5 (see above) 

MUSER 

ORDER 

PGSET RSTORE AUTOFL REQFLt 

PNCH 

PZSTRE 

t DYLOFLEX library routine 
CDC FORTRAN library routine 
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QCKROO 


Table 1. — (concluded) 
( AUTOFL REQFLt 


CLEAR 

MULLP 

. SECOND* 

QRROOT (see figure 19) 
READEM RDMAT 

READPT AUTOFL 

RLGRID / RWE12 

I RWI5 

RSTORE AUTOFL 


CCOMPE* 

CLDIV* 

CREV* 


READTPt 
TAPDIM 
^ WRITER 
REQFLt 
RUW (see above) 
RUW (see above) 
REQFLt 


RUW (see above) 


RWA10 

RWE12 

RWI5 

SENSOR 


TRPLOT 

USERS 

VREAD 


RUW (see above) 


RUW (see above) 
/ READTP* 
RWE12 
RWI5 

TAPDIM i 
AUTOFL 


RUW (see above) 
RUW (see above) 
ENDQR 
WARN 
REQFLt 


CLEAR 

RUW (see above) 


VUSER 

WA10 

WRITUM WRITER 

ZFIND RSTORE AUTOFL 

t DYLOFLEX library routine 
* CDC FORTRAN library routine 


AUTOFL 


REQFLt 


REQFLt 



2.1.3 OVERLAY QRM,7,0 


This overlay is the last logical step executed in the QR program. The maximum field 
length used during the QR execution is calculated and printed on the output file 
(TAPE6) and the dayfile, the SC4020 plot buffer is filled and emptied onto disk unit 99, 
local scratch files TAPE90, TAPE91 and TAPE96 are returned, and the CRU cost of the 
QR execution is calculated. 


Table 3. — Subroutines Called by Overlay (QRMJ,0) 


QRM70 calls 


BUFOUT 

' COST* 

EXIT* 

CLOCK* 

REMARK* 

TTY* 

QRCOST calls ^ 

PFAPPND* 

WA10 

GTACCID* 

RETURNFt 

FETADDt 

BUFFIL 

RETURNFt 


USERNUM* 


t DYLOFLEX library routine 
* CDC system library routine 


2.1.4 OVERLAY QRT,2,0 

This overlay calculates time response given a polynomial transfer function 
(ZEROS(s)/POLES(s)), time intervals, and print/plot instructions. All of this input 
information is read from data cards or calculated in the QRM0,0 overlay and passed 
into QRT,2,0 by named common. 


Sections 4 .3and 4.8 of volume I discuss the analytical approaches used to solve time 
response options available in QR. 

Once the time response has been calculated and stored in blank common, it is 
printed/plotted using the specifications established by the user. Section6.3.5of volume I 
details time response data cards. The time delay for the COMPUTE TIME RESPONSE 
WITH TAU is read and written from input and output respectively in this overlay. 

Figure 5 is a macro flow chart for the QRT,2,0 overlay. 

Table 4 displays a map of subroutines called in QRT,2,0. 
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Figure 5. - Macro Flow Chart of Overlay (QRT, 2,0) 




Table 4. — Subroutines Called by Overlay (QRT,2,0) 


OVERLAY (QRT,2,0) 
PROGRAM QRT20 


INLAP 

NEGMAT 

PLOTTM 


( 


/ 


CGLESMt 

GLESOMt 

CLEAR 

KGRID 


PLTLI 


PTPLOT 


SCALE1 

VPRNTD 

WARN 
WRTETPt 
/ PRPLOT 

^ SCALE2 


( 


RWE12 

TBLUn 

TPRINT 


DXDYV* 

GRID1V* 

SETMIV* 

IXV 

lYV 

LINEV 

A10R1 

BUFFIL 


BLKFLL 

PRLING 


t DYLOFLEX library routine 
CDC system library routine 



2.1.5 OVERLAY QRF,3,0 


This overlay calculates frequency and PSD response given a transfer function 
(ZEROS(s)/POLES(s)), gust spectrum, and print/plot instructions. All of this input 
information is read from data cards or calculated in the QRM,0,0 overlay and passed 
into QRT,2,0 by named common. 

Sections 4,6 and 4.7 of volume I discuss the analytical approaches used to solve for 
frequency and PSD response in the QR program. 

Once the frequency or PSD response has been calculated and stored in blank common, it 
is printed/plotted using the specifications established by the user. Section 6 . 3 • 6 of 
volume I details frequency response data cards. 

Figure 6 is a macro flow chart for the QRF,3,0 overlay. Table 5 displays a map of 
subroutines called in QRF,3,0. 



Figure 6. — Macro Flow Chart of Overlay (QRF,3,0) 
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Table 5. - Subroutines Called by Overlay (OR F, 3,0) 

Overlay (QRF,3,0) 

Program QRF3Q 


BDNYQ 

CLEAR 

FPRINT 


(see tree on following page) 


( CLEAR 
KGRID 


KPLOT 

PLTLT 


PRINTV 


^ DXDYV 
GRIDIV 
^ SETMIV 
POINTV 

( IXV 
lYV 
LINEV 


GMARGN 

MNSQVS 

PEAKS 

PMARGN 

RESP 


RW15 

SCALE2 

VPRNTD 

WARN 
^ WRTETPt 
RESP 


A10R1 

BUFFIL 


(same list as FPRINT above) 

RESP 

RESP 


t DYLOFLEX library routine 



Table 5. — (continued) 


BDNYQ Calls 

ARCPLT / PLTLI 


POINTV 
POLAR 
, POLER 


BRITEV 

CIRCIV / C1RC1V 


BRITEV 

SXDYV 

FAINTV 

FRAMEV 

LABARC 


RADI1V 


I XSCALV 
V YSCALV 


^ IXV 
lYV 

. LINEV 


^ PLTLI IXV 

lYV 

XSCALV 
, YSCALV 


^ LABLV BEDEXP 
PRINTV 

NXV 
, NYV 
^ IXV 
lYV 

LABLV (see above) 

NXV 

NYV 

^ PLTLI (see above) 


t DYLOFLEX library routine 



Table 5. ~ (concluded) 


BDNYQ 

CLEAR 

DRWGRD 

FAINTV 

KGRID 

POINTV 

PRINTV 

SETMIV 

VPRNTD 

WARN 

WRTETPt 


Calls 


LABLXC 

NODEC 

PLTLI (see above) 


( DXDYV 
GRID1V 
SETMIV 


A10R1 

BUFFIL 


t DYLOFLEX library routine 



2.1.6 OVERLAY QRR,4,0 


This overlay generates root locus plots given ZEROS(s), POLES(s), root locus 
eigenvalues, and plot specifications. All of this input information is read from data 
cards or calculated in the QRM,0,0 overlay and passed into QRR,4,0 by named common. 

Root locus eigenvalues are stored on the disk file TAPE96 and are read into blank 
common as soon as the overlay is entered. Once all of the root locus points have been 
read, each of the complex arrays corresponding to ZEROS(s), POLES(s), and root locus 
eigenvalues are sorted into descending order using the imaginary part of the arrays as 
the sort criteria. The arrays are then plotted in accordance to specifications established 
by the user. 

If QR is performing a COMPUTE FLUTTER AND ROOT LOCUS analysis during 
execution of this overlay, the ZEROS(s) and POLES(s) complex arrays are written to 
disk file TAPE96 for eventual processing in overlay QRR,5,0. 

Figure 7 is a macro flow chart for the QRR,4,0 overlay. 



Figure 7.— Macro Flow Chart of Overlay (QRR,4,0) 


Table 6 displays a map of subroutines called in QRR,4,0. 

2.1.7 OVERLAY QRR,5,0 

This overlay generates V-g plots given eigenvalues, k-values, reduced frequency 
reference length, and plot specifications. All of this input information is read from data 
cards or calculated in the QRM,0,0 overlay and passed into QRR,5,0 by named common. 




Table 6. — Subroutines Called by Overlay (OR R, 4,0) 

Overlay (OR RAO) 

Program QRR40 


CLEAR 

FSFt 

KGRID 


KPLOT 

PRINTV 

PRPLOT 

SORT 

VPRNTD 


^ DXDYV 
GRID1V 
. SETMIV 
POINTV 


( 

( 


BLKFLL 

PRLINE 

A10R1 

BUFFIL 


WARN 

WRTETPt 

t DYLOFLEX library routine 


Eigenvalues are stored on the disk file TAPE96 and are read into blank common as soon 
as the overlay is entered. As each set of eigenvalues is read, arrays corresponding to 
velocity (knots or kilometers) and damping are calculated. Once all of the eigenvalues 
have been read, the arrays of velocity, damping, and cycles per second are sorted into 
ascending order using velocity as the sort criteria. The arrays are then plotted in 
accordance to specifications established by the user. 

Figure 8 is a macro flow chart for the QRR,5,0 overlay. 



Figure 8. — Macro Flow Chart of Overlay {OR R, 5,0) 
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Table 7 displays a map of subroutines called in QRR,5,0. 


Table 7. — Subroutines Called by Overlay (OR R, 5,0) 


Overlay (QRR, 5,0) 
Program QRR50 


CLEAR 

FSFt 

KGRID 


^ DXDYV 
GRID1V 
< SETMIV 


PRINTV 

PRPLOT / BLKFLL 

I PRLINE 

STRMOVt 

VPRNTD / A10R1 

I BUFFIL 

WARN 

WRTETPt 


t DYLOFLEX library routine 


2.1.8 OVERLAY QRU,5,0 

This overlay (see table 8) reads data associated with the UNSTEADY AIRFORCE 
COEFFICIENTS data card. Section 4.9 of volume I explains in detail the unsteady 
aerodynamic options in QR. 


Table 8. - Subroutines Called by Overlay (QRU,5,0) 

Overlay (ORU,5,0) 

Program ORU50 

RWE12 

RW15 
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2.1.9 OVERLAY QRU,6,0 


The purpose of this overlay (see table 9) is to set up calls to overlays QRU,6,1; QRU,6,2; 
QRU,6,3; and QRU,6,6. 

Table 9. — Subroutines Called by Overlay (QRU,6,0) 

Overlay (QRU.6,0) 

Program QRU60 

AUTOFL REQFLt 

FSRt 

t DYLOFLEX library routine 


2.1.10 OVERLAY QRU,6,1 

The purpose of this overlay (see table 10) is to read data cards associated with the 
COMPUTE UNSTEADY FREQUENCY RESPONSE data card and increment the 
various loop counters associated with the unsteady frequency response option. Section 
if . 9 of volume I explains the detail of this QR option. 

Table 10. — Subroutines Called by Overlay (QRU,6,1) 

Overlay (QRU,6,U 
Program QRU61 

ATMS62 - 

ENDQR 

RWI5 


2.1.11 OVERLAY QRU,6,2 

The purpose of this overlay is to calculate the lower and upper frequency limits of the 
unsteady frequency response and to establish an array of k- values that can be used 
during the unsteady frequency response. 

Using the first aerodynamic influence coefficient matrix, its k-value, and equation 
if.9"l from volume I, a set of eigenvalues can be calculated that identify the range of 
frequencies associated with an aircraft. The eigenvalue array is scanned for the 
minimum and maximum non-zero roots. The lower and upper frequency bounds are set 
one decade lower and higher from the minimum and maximum complex roots. 

Once the frequency limits have been calculated, QR reads the specified disk containing 
aerodynamic influence coefficient matrices for the k-value associated with each of these 
matrices. The k -values are stored in an array for use in overlay QRU,6,3. 
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Figure 9 is a macro flow chart for the QRU,6,2 overlay. 



Figure 9.— Macro Flow Chart of Overfay(QRU,6,2) 


Table 11 displays a map of subroutines called in QRU,6,2. 

2.1.12 OVERLAY QRU,6,3 

The purpose of this overlay is to calculate the unsteady frequency response of equation 
14 . 9-1 in volume I. The user must provide the force vector array F and the normal 
coordinate position q as well as a number of other related inputs. Section k ^3 of 
volume I explains the detail of these inputs. Much of the input is read from data cards 
and stored in named common before this overlay is called. 

The frequency response itself is calculated by substituting j into equation I 4 . 9 -I and 
solving the resultant set of complex simultaneous equations for the user-identified 
normal coordinate position. Additional values of omega are picked in an ascending order 
such that the db change and phase change are less than the user-specified limit. As 
omega ascends in value, it is necessary to calculate the k-value corresponding to it. If 
the calculated k-value is greater than the k-values associated with the two aerodynamic 
influence coefficient matrices residing in core, new matrices are read from disk until an 
interpolation for the k-value calculated from omega can be made. 

The user may use disk numbers 1,2,3,4,8,9,11-20 to input the aerodynamic influence 
coefficient matrices. 


22 






Table 11.— Subroutine Called by Overlay (QRU,6,2) 


Overlay (QRU,6,2) 
Program QRU62 


AUTOFL 

REQFLt 


BSF* 

EPLUSD 

MOMEGA 

ORDER 

OR ROOT (see following page) 


RDMAT ^ 

[ READTP 



TAPDIM 

ENDQR 

RPRINT 

[ SECOND* 

WARN 

RSTORE 

1 STRMOVt 
AUTOFL 

RECFLt 

WARN 

t DYLOFLEX library routine 


* CDC FORTRAN library routine 




Table 11. — (continued) 


Routines called by QRROQT 
AUTOFL REQFLt 

CHESS / CLEAR 
I WRITUM 

CQREIG 


CSAVE2 

EIGMAT / EXPAND 


HELPER 


/ 


NORMAL 


/ 


GAUSPM AUTOFL 

PRECON / EQVCOL 
I IRSCAL 

RHESS 

RQREIG 

t DYLOFLEX library routine 
* CDC FORTRAN library routine 


WRITER 


CLEAR 

ROTSM 

WRITUM 

CLEAR 

COFSIN 

COLCHK 

ELMCHK 

POLCHK 

REGRUP 

ROOTS 

WRITUM 

CLEAR 

CLIM 

FINDPV 

GAUSS 

JJFIND 

NEGMAT 

REGRUP 

SWAP 

WRITUM 

REQFLt 


WRITER 


CLIM 

MULLP* 

WRITER 


WRITER 



Table 11. — (concluded) 


RST0R2 

SECOND* 

SVFORM 


WRITUM 


AUTOFL REQFLt 


^ EXPAND (see above) 
NORMAL (see above) 
MATPRT WA10 

SETORD / AUTOFL 
I JJFIND 
^ SWAP 

WRITER 


REQFLt 


t DYLOFLEX library routine 
* CDC FORTRAN library routine 



Figure 10 is a macro flow chart for the QRU,6,3 overlay. 



Table 12 displays a map of subroutines called in QRU,6,3. 


Table 12, — Subroutines Called by Overlay (QRU,6,3) 


Overlay (QRU,6,3) 
Program QRU63 


AUTOFL 

lORESP 


REQFLt 

BSF* 

CGLESMt 
CLEAR 
ETIMEC 
^ RDMAT 


SUBJOM 

MOVEFR 


RDMAT (see above) 


RSTORE 

TRANS 


( 


AUTOFL 

CLEAR 

ENDQR 


READTPt 

TAPDIM / ENDQR 
[^WARN 


REQFLt 


t DYLOFLEX library routine 
CDC FORTRAN library routine 


1-20 
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2.1.13 OVERLAY QRU,6,6 


The purpose of this overlay is to calculate the maximums and minimums, gain, and 
phase margin and to print/plot the frequency response calculated in overlay QRU,6,3. 
Input data for this overlay are passed by named common from the QRM.O.O overlay. The 
frequency response is passed from the QRU,6,3 overlay in blank common. 

Figure 11 is a macro flow chart for the QRU,6,6 overlay. 



Table 13 displays a map of subroutines called in QRU,6,6. 

Table 13. — Subroutines Called by Overlay (QRU,6,6) 


Overlay (QRU,6,6) 

Program QRU66 

BDNYQ (see table 5) 

UGMARG TBLUIt 

UNSQVS (see calls to MNSQVS in table 5) 

UPEAKS 

UPMARG TBLUIt 

t DYLOFLEX library routine 

2.1.14 OVERLAY QRU,7,0 

The purpose of this overlay is to set up calls to overlays QRU,7,1; QRU 7 3- and 
QRU,7,5. 
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2.1.15 OVERLAY QRU,7,1 


The purpose of this overlay is to read data cards associated with the COMPUTE 
FLUTTER 2 data card, increment various loop counters associated with the option, and 
read the disk containing aerodynamic influence coefficient matrices for the k-value 
associated with each of these matrices. The k-value is stored in an array for use in 
overlay QRU,7,3. 

Figure 12 is a macro flow chart for the QRU,7,1 overlay. 



2.1.16 OVERLAY QRU,7,3 

The purpose of this overlay is to calculate the eigenvalues of the matrix polynomial 
given by equation 4,g-lin volume I. Eigenvalues are calculated for each k-value that 
is input by the user. 

Overlay QRU,7,1 establishes a table of k-values that are used for interpolation purposes 
in this overlay. QR positions the disk file containing the aerodynamic influence 
coefficient matrices and reads two of the matrices into central memory. The first user 
specified k-value is inspected and if it lies between the first two k-values that reside on 
disk, the eigenvalues are calculated. Otherwise, another matrix is read from disk. This 
procedure is repeated until eigied a k-value that is below or above the range of the table 
read in overlay QRU,7,1; a diagnostic is written onto the output file. In addition, the 
user must input k-values for flutter analysis in ascending order. 

Figure 13 is a macro flow chart for the QRU,7,3 overlay. 

Table 14 displays a map of subroutines calls in QRU,7,3. 

2.1.17 OVERLAY QRU,7,5 

The purpose of this overlay is to plot the V-g analysis produced by overlay QRU,7,3. 
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Table 14. — Subroutines Called by Overlay (QRUJ,3) 


Overlay {QRU,7,3I 
Program QRU73 


AUTOFL 

FLUTEI 


REQFLt 
/ AUTOFL 
BSF* 
ENDQR 
FORMPM 
ORDER 
OR ROOT 
RDMAT 

RPRINT 

RSTORE 
^ WRITUM 


REQFLt 


(see table 11) 


( 

( 


RPADTPt 

TAPDIM 

SECOND* 

STRMOVt 

AUTOFL 

WRITER 


( ENDQR 
WARN 


REQFLt 


t DYLOFLEX library routine 
* CDC FORTRAN library routine 
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Eigenvalues are read into blank common from disk file TAPE96 as soon as the overlay 
is entered. As each set of eigenvalues is read, arrays corresponding to velocity (knots or 
kilometers) and damping are calculated. Once all of the eigenvalues have been read, the 
arrays of velocity, damping, and cycles per second are sorted into ascending order using 
velocity as the sort criteria. The arrays are then plotted in accordance to specifications 
established by the user. 

Figure 14 is a macro flow chart for the QRU,7,5 overlay. 



Figure 14. — Macro Flow Chart of Overlay (QRUJ,5) 
Table 15 displays a map of subroutine calls in QRU,7,5. 


2.2 DATA BASES 

The data bases for QR include several input, output, and scratch files as well as a large 
number of FORTRAN common blocks. 

2.2.1 INPUT DATA 

Data is input to the program in three forms; data cards, magnetic files, and 
user-provided FORTRAN subroutines. 

2. 2. 1.1 Card Input Data 

For a complete description of card input data, see section 6. 3 of QR usage document, 
volume I. 




Table 15. — Subroutines Called by Overlay (QRU,7,5) 


Overlay (QRU,7,5) 
Program QRU75 


CLEAR 

KGRID 


PRINTV 

PRPLOT 

STRMOVt 

VPRNTD 


^ DXDYV 
GRID1V 
SETMIV 


( 


BLKFIL 

PRLINE 

A10R1 

BUFFIL 


WARN 

WRTETPt 

t DYLOFLEX library routine 


2.2. 1.2 Magnetic File Input Data 

Magnetic tape can be used to supply matrices required by QR. These files should be 
copied to disk before being used with the QR program. Permissible disk file names are 
TAPEl, TAPE2, TAPE3, TAPE4, TAPES, TAPE9, TAPEll, TAPE12, TAPE13, 
TAPE14, TAPE15, TAPE16, TAPE17, TAPE18, TAPE19, TAPE20. Sections6. 1 and6,3 
of the QR User’s Manual provide additional information about tape/disk input. 

2.2. 1.3 FORTRAN Input 

The elements of the square polynomial matrix and column vector array used by QR can 
be input through the use of user-supplied FORTRAN subroutines. Sec. 6.8 of the QR 
User’s Manual describes this input option. 

2.2.2 OUTPUT DATA 

The results of a program execution can be printed or written on magnetic tape. In 
addition, a plot tape (TAPE99), plot vector tape, and punched matrix cards can be 
generated by the program. 
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2.2.2. 1 Printed Output Data 

A short sample of QR printed output is included in section 4.9 of volume I. 

2.2 .2 .2 Punched Card Output Data 

The PUNCH option described in section 6. 3. 2of volume I allows the user to punch the 
square polynomial matrix used by QR with the format 4(312, E 14.7). The matrix punched 
on cards may be read later by QR using the MATRIX or REPLACE card input options. 

2.2.2.S Magnetic File Output Data 

Two magnetic disk files are often produced by QR executions. TAPE92 contains data 
records that represent data appearing on SC4020 plots. 


TAPE99 is a binary file of plot instructions and data points for the SC4020 plotter. 
Users normally pass this information to the SC4020 plotter through use of the PLOTFIL 
command. Reference 3 provides additional information on this command. 

2.2.3 INTERNAL DATA 

Two methods are used to pass data from one portion of the program to another - 
common blocks and temporary disk storage files. 

2.2.3. 1 Common Blocks 

LABELED common blocks are used for communication between subroutines and 
overlays. The block names and contents are described on the following pages. Table 16 
shows the overlays in which each block is defined. 

BLANK common is used in secondary overlays as a variable length working area. 
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QRU.7,3 







Table 16. — (concluded) 


COMMON 

BLOCK 

NAME 

0'0'lAldD 

cc 

a 

cc 

a 

o 

CN 

H 

CC 

a 

o 

CO 

u." 

oc 

O 

CC 

CC 

a 

o 

tn 

cc 

cc 

O 

o 

in 

D 

CC 

a 

co“ 

5 

cc 

O 

co' 

D 

OC 

O 

CN 

(£ 

D 

OC 

O 

n 

to* 

D 

OC 

O 

CD 

cd' 

D 

OC 

a 

NUSFRP 









X 

X 


X 

X 

OMEGAS 









X 


X 

X 


OOS 












X 


OUTl 

X 



X 

X 

X 

X 






X 

PHASE 

X 













PLOT 

X 





X 

X 







PLOTFL 

X 






X 







PLOTFR 

X 




X 








X 

POLE 

X 





X 








PPLOY 

X 













PSDPLT 

X 




X 








X 

PTITLE 

X 



X 

X 

X 

X 






X 

PTRPLT 












1 


RCM 

1 

' X 











X 


ROWCOL 

1 X 













RWBLCK 

X 













RWBUFF 

X 

X 


X 

X 

X 

X 




X 


X 

SCLXX 

X 




X 








X 

SCLYY 

X 




X 








X 

TELEX 

X 

X 


X 

X 

X 

X 



X 




TIMES 

X 



X 










TIMOLD 

X 













TITL 

X 



X 

X 

X 

X 






X 

TPJUNK 

X 



X 





X 

X 

X 

X 

X 

TPJNK 





X 

X ; 

X 







UCONST 









X 

X 

X 

X 


UNITS 

X 






X 

X 


X 




UNSTDY 

X 







X 

X 

X 

X 

X 


UPOtNT 









X 

X 

X 

X 

X 

USERK 














VECTOR 

X 











X 


VTS 

X 



X 


X 

X 







ZERO 

X 





X 








ZPOLY 

X 













BLANK 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


X 


X 


X 


X 


X 

X 


X 


X 


X 


X 



X 



X 


X 


X 


X 


X 


X 


X 


QRU.7,0 




LABELED COMMON NAME: ALPHT 


DESCRIPTION; Tabic ofSC4020 hollerith characters 



* Refers to variable type: I = Integer 
R = Real 
C = Complex 
L = Logical 
H = Hollerith 

















LABELED COMMON NAME: CPOLl: 


DESCRIPTION: Array of noncancclled poles. 


Number 


Variable 


Dimension 


Engineering 

Nomenclalure 


Description 


1 NCPOLE 


Number of noncancclletl poles 


CPOLE 


C 


()0 


Eigenvalues that are desienalcil 
POLES. 






















LABELED COMMON NAME: DCGAIN 


DESCRIPTION: DC cdiii of feedback control system, etc. 


Number 

Variable 

T* 

Dimension 

Engineering 

Nomenclature 

Description 

1 

DCGAIN 

R 

1 


DC gain of feedback control system 


NINTEG 

1 

I 


Number of poles at origin user wishes 
to add to feedback control system 

3 

NDIFFR 

I 

1 


Number of zeros at ongin user wishes 
to add to feedback control system 

4 

ISTORE 

1 

1 


Variable flag QR uses to designate 
POLES. ZEROS, ROOT LOCUS, etc. 


* Refers to variable type: I = Integer 

R = Real 
C = Complex 
L = Logical 
H = Hollerith 




LABELED COMMON NAME: DEF 


DESCRIPTION: CanccllJtion tolerances 



Refers to variable Ivpe: I = Integer 

R = Reaf 
C = Complex 
L = Locical 
H = Hollerith 
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LABELED COMMON NA.ME; DUMBS 


DESCRIPTION; An-jy of SC4020 instructions 


Variable T* Dimension Descript 

Nomencbture ' 




LABELED COMMON NAME: DYN 

DESCRIPTION: Block of matrix dimensions and pointers 


Number 

Variable 

T* 

Dimension 

Engineering 

Nomenclature 

Description 

1 

MATORD 

1 

1 


Number of rows in A matrix 


MATDFl 

1 

1 


Highest power of S plus 1 

3 

NGRCOL 

1 

1 


Start of COL array in RWBUFF 

4 

NGRGRW 

I 

1 


Stan of GROW array in RWBUFF 

5 

NGRGCL 

I 

1 


Start of GCOL array in RWBUGG 

6 

NC.RIPV 

1 

1 


Start of IPIV array in RWBUFF 

7 

NGRNGR 

1 

1 


Two times number of guessed roots 

8 

NGRID 

1 

1 


Start of ID array in RWBUFF 

9 

NGROOT 

1 



Start of ROOT array in RWBUFF 


* Refers to variable type; I = Integer 
R = Reaf 
C = Complex 
L = Logical 
H = Hollerith 




LABELED COMMON NAME; FCONST 

DESCRIPTION: Block of constants for COMPLETE FLUTTEI 

^ : 

Number 

Variable 

T* 

Dimension 

Nomenclature 

Description 

1 

lALT 


1 

Counter of altitude positions 

-> 

NK 

I 

1 ! 

Number of k -values 

3 

NOEFK 

1 

1 

Number of k-value matrices to read 

4 

MN 

1 

1 

k -value matrices 


Refers to variable type; 1 = Intecer 
R = Real' 

C = Complex 
L = Logical 
H = Hollerith 




















LABELED COMMON NAME: FMTLBL 


DESCRIPTION: Formal orSC4020 label values 


Nuinber 

Variable 

■I 

Dimension 

Engineering 

Nomenclature 

Description 

1 

IFMT(2) 

I 



Formal of SC4020 label value 


LMARY 

I 

’ 



3 

LADJX 

1 

1 


Constants for label positioning. Sec 
subroutine NODEC 

4 

LBMARX 

I 

1 




* Refers to variable tvpe: I = Inteeer 

R = Real' 

C = Complex 
L = Lot 

H = Hoi 
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LABELED COMMON NAME: FRMN 


DESCRIPTION: Frame count 



Refers to variable type: I = Integer 
R = Real 
C = Complex 
L = Locical 
H = Hollerith 









♦ Refers to variable type: I = Inlecer 
R = Reaf 
C = Complex 
L = Logical 
H = Hollerith 


54 















LABELED COMMON NAME: IGARFR 


DESCRIPTION: Block of varijhlc 

‘s used in ur 

steady aerodynamics options 

Number 

Variable 

D 

Dimension 

Engineering 

Nomenclature 

Description 


GINVRS 

R 

6 

i/7 

Inverse gamma values 

■» 

ALTITU 

R 

6 


Altitude positions for flutter 
analysis 

3 

RFRL 

R 

I 

b 

Reduced frequency reference 
length 

4 

NIGFAA 

I 

1 


Number of inverse gamma values 

5 

NIG 

I 

1 


Current inverse gamma value 
subscript 

6 

NKVALU 

I 

I 


Number of k -values on disk 

7 

VALUER 

R 

1 


Current k-value 

8 

1 

NK 

I 

1 


k*value counter 

1 















LABELED COMMON NAME. ILOW 


DESCRIPTION; Hollerilh table of OR commands 



R = Re; 
C = Coi 
L = Loi 
H = Hoi 
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LABELED COMMON NAME; ITEM 


DESCRIPTION: CRU cost tracking information 



* Refers to variable type: 1 = Integer 
R = Reaf 
C = Complex 
L = Logical 
H = Hollerith 






LABELED COMMON NAME; JKL 


DESCRIPTION; Lcydinii coiTficients of POLE and ZERO polynomials 


Number 

Variable 

T* 

Dimension 

Engineering 

Nomenclature 

Description 

1 

PCOEF 

C 

1 


Leading coefficienl of POLE 
polynomial 


2COEF 

C 

I 


Leading coefficienl of ZERO 
polynomial 











LABELED COMMON NAME; MATCOM 


DESCRIPTION: Constants used- for ciiicnvalue rootin:^ 


Number 

Variable 

T* 

Dimension 

Engineering 

Nomenclature 

Description 


COFF 

C 

1 


Leading coefficient of A matrix 


PIVOT 

C 

1 


Value of last pivoted point referred to 

3 

MORIG 

I 

1 


Original number of rows in A 
matrix 

4 

NORM 

L 

1 


NORM is TRUE, if subroutine 
CLIM is to do Gaussian elimination. 

5 

SINGLR 

L 

1 


SINGLR is .TRUE, if A is sigular 

6 

VTEST 

L 

1 


VTEST is .TRUE, if state variable 
expansion is underway 

7 

ZTEST 

R 

1 


A number of considered to be 
zero if less than or equal ZTEST 

8 

EPS 

R 

1 


Two polynomials are equal if their 
normalized coefficients vary by no 
more then EPS 


* Refers to variable type: I = Integer 
R = Real 
C = Complex 
L = Logical 
H = Hollerith 










LABEL! 

DESCRI 

ID COMMON NAM 
PTION; Muximu 

E: MA 

n frequ 

XFR 

ency of root 

s printed by OR 


Number 

Variable 

T* 

Dimension 

Engineering 

Nomenclature 

Description 

1 

MAXFR 

R 

1 


Eigenvalues having a complex 
absolute value greater than 
MAXFR will not be printed by QR 
during a KIT execution 


* Refers lo variable type: 1 = Inieeer 
R = Reaf 
C = Complex 
L = Logical 
H = Hollerith 
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LABELED COMMON NAME: 


MSZ 


DESCRIPTION Si/e ol mutrix stored on TAPE 90 and TAPE 91. 



* Refers to variable type; I = Integer 
R = Real 
C = Complev 
L = Logical 
H = Hollerith 


3 
















LABELED COMMON NAME: NFF 


DESCRIPTION: Number of forcing function points 



* Refers to variable type; 1 = Inteeer 
R = Reaf 
C = Complex 
L = Locical 
H = HoUerith 


5 










R = Real 
C = Complex 
L = Logical 
H = Hollerith 

























LABELED COMMON NAME: PLOT 


DESCRIPTION: Root locus plot specifications 


Number 

Variable 

T* 

Dimension 

Engineering 

Nomenclature 

Description 

1 

NPLOT 

1 

1 


Number of root locus regions 

•> 

XLEFT 

R 

6 


Minimum values of real axis 

3 

XRT 

R 

6 


Maximum values of real axis 

4 

YTOP 

R 

6 


Maximum values of imaginary axis 

5 

YBOT 

R 

6 


Minimum value of imaginary axis 

6 

NPHA 

I 

1 


Flag indicating if SC4020 plots 
phase locus are to be separate 

7 

IRLPF 

I 



Flag indicating if SC4020 plots 
are to be produced. 

8 

NRNTS 

I 



Number of loci sets on TAPE96 

9 

MFILE 

1 



Number of loci files on TAPE96 

10 

NRNTSS 

I 



Count of loci sets on TAPE96 
for flutter locus 


* Refers to variable type: I = Inlefier 
R = Reaf 
C = Complex 
L = Logical 
H = Hollerith 









LABELED COMMON NAME: PLOTFR 

DESCRIPTION: Plot specification for Bode and Nyquist Plots 


Number 


Variable 

NPLOTF 

IFROPT 

lOMGAL 

lOMGAR 

IBDNYO 

IDBSCL 

ILEFTW 

IRIGTW 


Dimension 


Engineering 

Nomenclature 


I 

1 

b 

6 

6 

6 

6 

6 


Description 


Number of regions to plot 

Flag indicating automatic 
scaling 

Log to base 10 of left limit of 
frequency scale. 

Log to base 10 of right limit 
of frequency scale. 

Flag indicating Bode and/or 
Nyquist plots. 

Flag indicating automatic 
scaling 

Temporary location for left 
scale 

Temporary location for right 
scale 



















LABELED COMMON NAME: PSDPLT 


DESCRIPTION: PSD plot 

specifics 

lions 



Number 

Variable 

T* 

Dimension 

Ejigineering 

Nomenclature 

Description 

1 

4 

NPSPP 

ISCPS 

FLR 

lORDS 

I 

1 

R 

H 

I 

I 

i: 

8 


Flag indicating if there are to 
be PSD plots 

Flag indicating cps or radians 
along abscissa 

Maximums and minimiims of abscissa 
Optional title for PSD plot 
































LABELED COMMON NAME: RCM 


DESCRIPTION: Ckuii locjtion informalion 



* Refers to variable tvpe: 1 = Integer 
R = Reaf 
C = Complex 
L = Logical 
H = Hollerith 




LABELED COMMON NAME: ROWCOL 


DESCRIPTION; Logical flag that indicates rows or columns of the A matrix 


Number 


Variable 


T* 


Dimension 


Engineering 

Nomenclature 


Description 


I IRC 


L 


If IRC is .TRUE., scanning a row. 
.FALSE, means column. 


Refers to variable type: I = Integer 
R = Real 
C = Complex 
L = Locical 
H = Hollerith 
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* Refers to variable type. I = Integer 
R = Reaf 
C = Complex 
L = Logical 
H = Hollerith 
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LABELED COMMON NAME: TPJUNK and TPJNK 


DESCRIPTION; Scratch array for RLADTP and WRTE:TP subroutines 


Number 


Variable 


Dimension 


Engineering 

Nomenclature 


Description 


1 

3 

4 

5 

6 


NTAPE 

NFILE 

NMAT 

MATDC; 

NPRT 

NAME 


Unit number 
Number of files to skip 
Number of matrices to skin 
Power of S matrix plus 1 
Print flag 

Fixed point number designatin 


7 M II 

S N II 

9 B I i: 

10 IRROR I 1 


matrix 

Row dimension 
Column dimension 
Identification array 
Error return 





LABELED COMMON NAME: UNCONST 


DESCRIPTION: Unstcjdy aerodynamics information 



Refers to variable type: I = Integer 
R = Real 
C = Complex 
L = Logical 
H “ Hollerith 






LABELED COMMON NAME: UPOINT 


DESCRIPTION; Unslcjdy aerodynamics information 


Number 

Variable 

T* 

Dimension 

Engineering 

Nomenclature 

Description 

1 

— 

lALT 

I 

1 


Altitude array dimension 


NMATR 

1 

1 


Number of matrices to forward 
space on first k -value 

3 

NK 

‘ 

1 


Number of k -values 

4 

NCLOL 

1 

1 


Flag that determines kind of 
frequency response 

5 

NFRLIM 

j 

1 


Flag that determines detail 
of frequency response 

6 

NEOFK 

i 

i 

I 


Flag that determines if air 

force matrices have been normalized 

7 

NXYMI 

1 

1 

j 

I 

1 


Number of combinations for 
closed loop response analysis 

8 

IXYMI 

! 

1 


Loop counter for NXYMI 

9 

IRHSV 

.1 

I 

1 


Loop counter for number of 
vectors 


* Refers to variable type: I ” Integer 

R *• Real 
C - Complex 
I = Logical 
H - Hollerith 




LABELED COMMON NAME: UNSTDY 


DESCRIPTION; Unsteady aerodynamics information 


Number 

Variable 

T* 

Dimension 

Engineering 

Nomenclature 

Description 

I 

MACHNO 

R 

1 


Mach number 

-> 

RFRL 

R 

I 

b 

Reduced frequency reference 
length 

3 

CF 

R 

1 


Correction factor 

4 

RHO 

R 

1 


Density of air 

5 

V 

R 

1 


Free stream velocity 

6 

NALT 

I 

I 


Number of altitudes 

7 

ALT 

R 

6 


Altitude values 

8 

LRHOAV 

L 

I 


If LRHOAV is .TRUE., RHO and 
V are to be calculated bv 
ATMS62 


• Refers (o variable type: I Integer 
R = Reaf 
C = Complex 
L = Logical 
H = Hollerith 




LABELED COMMON NAME: USERK 


DESCRIPTION: L)scr-spi*cificd K-valucs 



Refers to variable type: I = Integer 

R = Real 
C = Complex 
L = Locical 
H = Hollerith 
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LABELED COMMON NAME: VTS 


DESCRIPTION; Onlinatc title inrornution for this response 



Refers to variable type: I = Integer 
R = Reaf 
C = Complex 
L = Logical 
H = Hollerith 






labeled common NAME: ZERO 


DESCRIPTION: Amy of user indii 

;ated ZEROS 

Number 

Variable 

T* 

Dimension 

Engineering 

Nomenclature 

Description 

I 

NZERO 

ZERO 

I 

L 

1 

60 


Number of eigenvalues in ZERO 
array 

Eigenvalues that have been 
designated ZERO(S) 




LABELED COMMON NAME; 


ZPOLY 


DESCRIPTION: Polynomial coefficients of ZERO(S) array 



* Refers to variable type: I = Integer 

R = Real 
C = Complex 
L = Logical 
H = Hollerith 
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2.2.S.2 Magnetic Files (Disk) 

QR uses three disk files for temporary storage of data. Data is written and read from 
two of these files using BUFFER OUT and BUFFER IN respectively. The third file is 
written and read using unformatted FORTRAN WRITE and READ operations. 

TAPE90 Scratch File 

TAPE90 is a scratch file written and read by the QRM,0,0 overlay. This file contains an 
exact copy of the A matrix created by the user with the MATRIX/REPLACE/READ 
MATRIX TAPE etc. input data cards. 

The file structure is: 

Matrix Size Matrix 


NROWxNCOLxNDEGx2 


where: 


NROW 

= the number of rows 

NCOL 

= the number of columns 

NDEG 

= the highest power of S plus 1 
designates complex elements 



TAPE91 Scratch File 

TAPE91 is a scratch file v/ritten and read by the QRM,0,0 overlay. This file contains an 
exact copy of the matrix for which eigenvalues are to be calculated. If the algorithm for 
eigenvalues fails because of a singular matrix, QR will reread TAPE91, set a trace flag, 
and reattempt the eigenvalue solution printing traceback information on the output file 
until the singularity condition reoccurs. 

The matrix stored on TAPE91 has the same dimensions as the matrix stored on 
TAPE90. (However, element values may be different.) 


no 


TAPE96 Scratch File 


Tape96 is a scratch file written by the QRM,0,0 and QRU,7,3 overlays and read by the 
QRR,4,0, QR,5,0; and QRU,7,5 overlays. In addition, after the QRR,4,0 overlay has read 
the intended set of root locus points, it writes the open loop zeros and poles (from the 
root locus calculation) and an end-of-file onto TAPE96. 


The file structure is: 



The number of eigenvalues (NRl) 
k-value associated with eigenvalues 
starting in words. 

Eigenvalues in complex format. Two 
times NRl words 

The number of eigenvalues (NR2) 

k-value associated with eigenvalues 
following 

Eigenvalues in complex format. 

Two times NR2 words 


mber of eigenvalues, k-value, eigenvalues - is repeated until 
all flutter eigenvalue rootings have been completed. If the program is performing a 
COMPUTE FLUTTER AND ROOT LOCUS operation, the open loop poles and zeros are 
written onto TAPE96 followed by an end-of-file once the gain and/or phase loci have 
been read by the root locus plotting overlay (QRR,4,0). The program then advances to 
the next k-valuc and writes a second set of records onto TAPE96 followed by another 
end-of-file. This process is repeated until all k-value rootings have been completed. 

In the case of a COMPUTE FLUTTER, COMPUTE FLUTTER 2 or COMPUTE ROOT 
LOCUS, only one file is written. 


2.3 SUBROUTINES 

Short descriptions of the purpose of each function and subroutine used in QR follow in 
alphabetical order. 

Subroutine AA (OMEGA, CPS, FMAG, DB, PHI) 

AA sets up subroutine calls to print and plot frequency and PSD response, finds 
maximums and minimums, calculates gain margin, and phase margin for the unsteady 
frequency response option. 
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Subroutine ACCT 


ACCT contains a number of entry points localizing calls to QR overlay segments. Entry 
points for ACCT are FRESP, FTPLOT, TRPLOT, ENDQR. Before an overlay is called, 
subroutine AUTOFL is called to set the core size for the overlay execution. Only the 
overlays that are called from more than one location in QR are centralized in this 
routine. The overlays called are: QRM,1,0; QRF,3,0; QRR,4,0; QRR,5,0; and QRM,7,0. 

Subroutine ARCPLT (RAD, THETA, NUM, RI, AL, INC, IND, ISYM) 

ARCPLT is a routine to connect points on a polar grid by linear lines in the polar space. 
It also may place a symbol at the point and may connect the points with linear lines in 
the cartesian space. 

INPUT: 

RAD An array of radii in radial units 
An array of angles in degrees. 

RAD and THETA have 1 to 1 correspondence. 

NUM The number of points to process. 

^1 no processing will take place. 

AL Angular otTset in degrees. The value in the 3 o’clock position on the polar 

grid. 

INC Line density. 

< 1 Straight lines will connect the points. 

ARCS will be constructed between the points using chains of straight 
lines of INC angular length. 

IND Place a symbol on the point. 

^0 Do not place a symbol. 

>0 Do place a symbol. 

ISYM The number of the symbol to use. 

OUTPUT: None. 

Function ARCTAN (FRR,FRI) 

ARCTAN calculates the phase angle in degrees associated with the complex number 
FRR -f j FRI. The angle returned by ARCTAN will be in the range -180 to + 180 
degrees. 

Subroutine ARCIV (R, DR, M, LA, KDO) 

ARCIV was designed to be a primary level routine in the creation of a polar grid. Care 
has been taken to insure uniform density lines and to provide for minimum execution 
time. 
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ARCIV draws the axes of equal radii for the polar grid. 

INPUT: 

R The radial dimension of the circle. 

DR The radial dimension increment to be used in forming the circles. 

M The index to be used with DR to form the circles. Every Mth DR will be 

formed. 

LA Specifies the major axis of the grid. 

= 1 The vertical axis is major. 

^ 1 The horizontal axis is major. 

KDO The nature of the display. 

= 0 Only outermost circle will be drawn. 

0 Every Mth DR will be drawn along with the circle at R. 

OUTPUT: None 


Subroutine ASSIST 

ASSIST reads namelist data from the INPUT file and writes each card image read to 
the OUTPUT file and disk unit 7 (TAPE7).This subroutine is described in more detail in 
section 6.8 of volume I. 

Subroutine ATMS62 (ALT, SONIC, PRESS, DENSI, TEMP, MACHNO) 
ATMS62 calculates atmospheric properties of the US standard atmosphere, 1962. 

Subroutine AUTOFL 

AUTOFL calculates the field length required for execution and requests a field length 
change through subroutine REQFL. In addition, the maximum field length and last 
word address are calculated. These variables will be written out if the TRACE flag is 
".TRUE.”. 


Subroutine AXISD (IX, lY, IIX, 0) or (IX, lY, IIY, 1) 

AXISD draws a horizontal or vertical line from point (IX, lY) in raster coordinates to 
(IIX,IY) or (IX,IIY). Note that (IIX-IX) ^ 64 or (IIY-IY) 5= 64. 

Subroutine AlORl (A,B,N) 

AlORl takes N hollerith characters stored in array A and creates N "Rl” words in 
array B. 
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Subroutine BDNYQ (OMEGA, DB, PHI, N) 


BDNYQ is the QR subroutine that creates Bode and Nyquist on the SC4020. 
Section 4.7.6 of Volume I provides additional information on the QR frequency response 
option. 


Subroutine BEGIN 

BEGIN initializes default values for QR. See Sec^ 6,3^1 of Volume I for additional 
details. 


Subroutine BIGGER 

BIGGER is used to extend dynamic data arrays in C^R. 


BLOCK DATA 

BLOCK DATA initializes all named common blocks in QR. 

Subroutine BRITEV 

BRITEV sets the SC4020 exposure to heavy. In addition, entry FAINTV is included in 
this subroutine to set the exposure to light. 

Subroutine BUFFIL 

BUFFIL takes the current set of SC4020 instructions and packs them into the TAPE99 
output buffer OUTI. When OUTI is filled to 408 CDC words, subroutine BUFOUT is 
called to empty the buffer and reset the buffer count pointer. 

Subroutine BUFOUT 

BUFOUT performs a BUFFER OUT to transfer labeled common block OUTI, array 
lOUT, to disk unit 99. Once the BUFFER OUT operation is completed, counters are set 
to reflect that the disk unit 99 buffer is empty. 

BUFOUT is called to empty the buffer and reset the buffer count pointer. 

Subroutine CANCEL 

CANCEL takes the user designated POLE(S) and ZERO(S) arrays and determines if 
pole-zero pairs can cancel. The subroutines prints cancelled pole-zero pairs, remaining 
zeros and polynomial coefficients. 
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Subroutine CCOMPE (A,NA,CX,CY) 

CCOMPE evaluates a polynomial with complex coefficients given in the array A. The 
complex value returned is CY for the complex point CX. 

Subroutine CHESS (F,MS,B) 

CHESS forms the upper hessenburg matrix of the input matrix F, which is MS by MS. 

The hessenburg form is stored over the input matrix F. B is a scratch complex vector. 

Subroutine CIRCIV (L,RI,RO,Al,Dl,KA,KODE,KDO,AMM) 

CIRClV’s purpose is to draw circular grids on the SC4020. 

CALL CIRCIV (L,RI,RO,Al,Dl,KA,KODE,KDO,AMM) 

INPUT: 

L Frame advance control. 

do not advance frame. 

^0 advance the frame. 

RI Value in radial units of the center of the displayed grid (radial offset). 

RO Value in radial units of the outside of the grid. 

A1 Angle to be displayed at the 3 o’clock position on the grid, modulo 360 (Angular offset). 

D1 In a counterclockwise direction the angle difference used for display of radial lines. 

D1 = 0 no lines drawn. 

KA The direction of plot angles. 

= 1 The angles appear in a counter clockwise direction. 

= 1 The angles appear in a clockwise direction. 

KOBE Grid type control. 

= 1 Circular. 

= 2 Left semicircle displayed at right hand edge. 

= 3 Right semicircle displayed at left hand edge. 

= 4 Top semicircle displayed at bottom edge. 

= 5 Bottom semicircle displayed at top edge. 

= -1 Circular with entire grid within margins. Allows truncation of the grid if 

the margins do not specify a square space. 

KDO Grid Control, 

= 0 only the outer most circle will be displayed with no labeling, 

^ 0 Full gradation will occur along with full labeling. 
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AMM Raster increment between minor arcs usually between 8. and 30 (See 
DXDYV). 

OUTPUT: None 


STATUS 
AFTER CALL: 

Scale factors will be set up for a polar mapping onto a cartesian grid, center of 
the plot of this grid is at point (0,0) and the outside of the circle is at a 
distance (RO-RI) from the center. In the case of a circular plot, the point (0,0) 
will be at the center of the space set up by the margins; in the case of a 
semicircle, it will be at the center of the margin at the open end. In all cases, 
the angle that will be mapped at the 3 o'clock position is 0® = (Al- Al).The 
space is set up so the angles when mapped will be ascending in the clockwise 
direction for clockwise space; descending in the clockwise direction for a 
counterclockwise space. All radii drawn will be labeled. To suppress radial 
lines, set D1 = 0. 


Subroutine CLDIV (A,N,B>C,M) 

CLDIV divides a complex polynomial by a linear expression (Z/B). A is array 
of complex coefficients; N is degree of polynomial; B is complex constant in 
linear factor (Z + B); C is complex array containing quotient and remainder. 
C(l) contains the remainder and C(2) through C(N + 2) contain the 
coefficients of the quotient polynomial contained in C(2) through C(M + 2). 

Subroutine CLEAR (M,N) 

CLEAR zeros N elements in array M. Entry point NEGMAT reverses the sign 
of N elements in array M and entry BLKFLL fills N elements in array M with 
blank characters. If N is less than or equal zero, no action is taken. 

Subroutine CUM (U,MATORD,MATDEl,II,JJ,KK,ROW,COL) 

CLIM performs gaussian elimination, by rows or columns. Entry point GAUSS assumes 
the pivotal element has been identified. 

Subroutine COFSIN (II, ROW) 

COFSIN recalulates the leading coefficient of the matrix polynomial based on the 
number of permutations of rows or columns. 

Subroutine COL 

COL calculates the branch flag ICOMTY. This branch flag is then used in the QR main 
program. 
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Subroutine COLROW (U,MATORD,MATDEl) 


COLROW is called when the user has requested a PRINT PARTIAL MATRIX input 
option (see sec. 6.3.2 of Vol . I). 

Subroutine CQREIG (A,N,ROOT,V,U) 

CQREIG computes the eigenvalues through use of the QR similarity transformation. A 
is an N by N matrix; ROOT is the eigenvalue array; and V and U are scratch vector 
arrays. 


Subroutine CRAMER (U, VECTOR, MATORD,MATDEl) 

CRAMER substitutes the VECTOR array into a user defined column of the matrix 
polynomial U. MATORD and MATDEl are dimension specifications: 

Subroutine CREAD (U,MO,MD,V,ME,IROW,ICOL) 

CREAD performs the input options of CONTINUOUS MATRIX, DELETE, DELETE 
AND REDUCE, MATRIX, MATRIX COMPLEX, REPLACE, REPLACE COMPLEX, 
TRUNCATE, and TRUNCATE and REDUCE. Further information on these program 
options can be obtained from section 6.3.2 of Volume I. 

Subroutine CREV (P,NP,PR) 

CREV reverses the order of coefficients of a complex polynomial in an array. P is the 
array of polynomial coefficients which will be reversed; N is degree of polynomial stored 
in P; and PR is array of polynomial coefficients stored in reverse order. 

Subroutine CSAVE (U) 

CSAVE saves the matrix polynomial on disk unit 90 after certain program options have 
been executed. For example, after data associated with the MATRIX input option has 
been read by QR, subroutine CSAVE is called. In addition, CSAVE contains the entry 
point CSAVE2. This section of the subroutine saves the matrix polynomial on disk unit 
91. CSAVE2 is called at the beginning of an eigenvalue rooting. 

Subroutine DXDYV (1,XL,XR,DX,N,I,NX,DC,IERR) 
or (2,YB,YT,DY,M,J,NY,DC,IERR) 

A call to this subroutine will help provide an optimum calling sequence for GRIDl V. Output 
from DXDYV will be DX,N,I,NX, and lERR or DY,M,J,NY, and lERR. Additional 
information on XL,XR,YB,YT,DX,CY,N,M,I,J,NX, and NY is located with Subroutine 
GRIDl V. DC limits the density of the grid. The grid lines drawn by GRID IV, using 
arguments supplied through DXDYV, will be no closer than DC raster positions. DC should 
never be less than 3.0. DC = 8.0 or 20.0 are recommended values. lERR is an error indicator. 
If it is non-zero, then subsequent use of GRIDl V will probably not produce a satisfactory grid 
with the given arguments. 
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Subroutine DRWGRD (IOMGAR,IOMGAL,DBTOP,DBBOT) 


DRWGRD creates the grid for bode frequency plots based on abscissa and ordinate 
minimums and maximums specified by the user. 

Subroutine EIGMAT (U, MATORD, MATDEl, ID, ROW, COL, GROW, GCOL, 

IPIV, ROOT) 

EIGMAT sets eigenvalue rooting constants and calls a number of subroutines to place 
matrix polynomial in normalized form. 

Subroutine ELMCHK (U,MATORD,MATDEl,II,JJ,KK,ROW,COL) 

ELMCHK scans the matrix for reducible rows or columns. U is the matrix array, MATORD 
is number of rows in U; MATDE 1 is highest power of S plus 1 . II, JJ and KK identify the row, 
column and S + 1 pivotal location for reduction. ROW and COL are logical arrays to keep 
track of which rows and columns have been scanned. ELMCHK is also entered at COLCHK 
and POLCHK. 


Subroutine EPLUSD (U,MATDURD,MATDE1) 

EPLUSD adds two matrices. It is used when calculating unsteady frequency response. 

Subroutine EQVCOL (N,IDIM,A,ISTACK,LOWST,LOC,IEQUIV,INIX) 

EQVCOL detects and isolates equivalence classed under reachability within a graph, 
given a connection matrix of the graph. 


Subroutine ETIMEC (AK2,VK) 

ETIMEC forms a new matrix from a straight-forward multiplication. 

Subroutine EXPAND (U,MATORD,F,DF,ISIZE,IORD) 

This subroutine expands a matrix polynomial that has second, third, etc. order terms 
(S , . , . ) to a first order system (S^ terms only). 

Subroutine FINDPV (U,MATORD,MATDEl,II,JJ,RDW,COL) 

FINDPV finds the complex element along a given row or column that has the largest 
magnitude. 


Subroutine FLUTEI (AK1,AK2,U,MX,NX) 

FLUTEI determines the two k-value matrices that must reside in core for the proper 
interpolation. The matrix polynomial corresponding to the user’s input k-value is then 
formulated, rooted, and printed. In addition, eigenvalues and k-values are placed on 
disk (TAPE96) for plotting. 
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Subroutine FORMPM (ROOT,NROOT,FUNCT,RFUNCT,COEF,NPOLY) 

FORMPM forms a polynomial array of NROOT eigenvalues, the array ROOT. COEF is 
the coefficient of the highest power of S. 

Subroutine FRAMEV (N) 

FRAMEV advances the SC4020 microfilm, places the data, page number and titles on 
the new grid. 


Subroutine FRGRID 

FRGRID reads input data associated with the PLOT FREQUENCY RESPONSE and 
GUST SPECTRUM commands. Section 6. 3* 6 of Volume I contains information about 
those data cards. FRGRID also has the entry point gust. 

Subroutine FRLINE (X.Y,M,SYMB) 

FRLINE is the same as PRLINE except that it is used in the COMPUTE FLUTTER 2 
overlay, QRU,7,5. 


Subroutine GAINS 

GAINS reads input data associated with the GAINS and DC GAIN commands. Section 
2.3.3. ,Vol. I, has additional information on these data cards. GAINS also has the 
entry point GAINDC. 

Subroutine GAUSPM (U,MATORD,MATDEl,MATONE) 

GAUSPM determines the highest power of S that is used in a matrix polynomial. 

Subroutine GRIDIV (L,XL,XR,YB,YT,DX,DY,N,M,I,J,NX,NY) 

This routine will produce a grid which has some lines emphasized and some lines 
labeled, as desired. (Margin spaces will be reserved at the top, left side, and bottom of 
the grid. These margins will normally be 24 raster counts wide.) The routine sets 
exposure to heavy. At the completion of the routine, scale factors will have been 
established and made available, for the conversion of floating point values to raster 
coordinates, during the operation of subsequent sub-programs in the package, 

L = 0 The film will be advanced, ID and corner marks will be printed. 

L = 1 The film will be advanced, the date and frame count will be placed in the 

upper right-hand corner of the frame. Corner marks will be suppressed. 

L = 2 The film will not be advanced and the identification information will not be 
displayed. 

L = 3 The film will be advanced but the identification information will be 
suppressed. Corner marks are printed. 
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L = 4 The film will be advanced and both the ID and the corner marks will be 
suppressed. 

XL,XR Floating point values of X for the left-most and right-most limits of the grid, 

YB,YT Floating point values of Y for the bottom and top limits of the grid. 

DX,DY Floating point data increments at which vertical and horizontal grid lines, 
respectively, will be displayed. If DX = 0, then no vertical grid lines will be 
shown. To prevent crowding, no more than 120 grid lines in one direction will 
be drawn. 

Integers that cause every Nth vertical grid line and every Mth horizontal grid 
line to be retraced for emphasis. If N = 0, then no vertical line is emphasized. 
If M = 0, then no horizontal line is emphasized. 

If either N or M is negative, a square grid will be drawn. Note that if N = 0 
or M = 0 and a square grid is required, then both N and M should be made 
negative. 

Integers that cause every Ith vertical grid line and every Jth horizontal grid 
line to be labeled. If I = 0 then no vertical line is labeled. If J == 0 then no 
horizontal line is labeled. Note that label margin space is in addition to the 
margin reserved for titles. 

Integers indicating number of characters to be displayed in the labels of 
vertical and horizontal lines respectively, 

A positive sign will give a label in decimal format. A decimal point must be 
counted as one of the characters, but the sign is not counted. The largest 
number of digits permitted is 6; i.e., NX(or NY) can be + 7 at rnost, A 
negative sign will produce a label in scientific notation (e.g., 1.25x10^), The 
sign, decimal point, and exponent will be displayed in addition to NX (or NY) 
characters. NX (or NY) must be no less than -6. 

Subroutine GMARGN (N) 

GMARGN calculates the gain margin of an input frequency response. 

Subroutine HELPER (U, MATORD,MATDEl,MATSIZ, ROW, COL, ADDRT, TEMP) 

HELPER will compress the input matrix polynomial by: removing rows or columns that 
are zero except for one element, removing rows or columns with like powers of S, 
removing rows or columns with a common polynomial factor. Additional information on 
these reduction techniques is available in reference 2 , 

Subroutine INLAP (NT) 

INLAP calculates the time response for an expression in Laplace transform with 
multiple poles. Immediate calculation by residue method is performed for single poles, 
while for the remaining poles, calculation is done by solving a linear system of 
equations. Further technical information on INLAP is presented in section 
Volume I. NT is the number of time history data points. 


N,M 


I,J 


NX, NY 
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Subroutine lORESP (U,CVEC,IPR,MATORD,MATDEl,AK2, SO, OMEGA, 

DB,PHI,VECTOR,CVECT) 

lORESP calculates the unsteady frequency response. 

Subroutine IRSCAL (A,T,IA,N,M) 

gIRSCAL provides scaling for irreducible submatrices. 

Function IXV(X) 

IXV converts a floating point abscissa value X to a fixed point raster value IX ready to 
plot. A previous call to XSCALV must have been made before using IX = IXV(X). IX = 
A*X = B. 


Function lYV(Y) 


lYV is similar to IXV, setting Y for X, lY for IX, C for A, D for B and YSCALV for 
XSCALV. 


Function JJFIND(IPIV,II,N) 

JJFIND finds the row/column corresponding to pivot located in column/row II. 

Subroutine KGRID (LADV,XL,YB,XR,YT,LOG) 

KGRID sets the raster positions for linear or logarithmic plots, calls subroutines 
DXDYV and GRIDIV based on the input LOG (LOG = -1, time response grid; LOG = 0, 
root locus grid; LOG = 1, frequency response grid). XL and XR are maximum and 
minimum abscissa values, YB and YT are maximum and minimum ordinate values. 

Subroutine KPLOT (X,Y,IPTS,KHAR,XLF,YBOT,XRT,YTOP) 

KPLOT plots the character specified in KHAR for the X-Y pairs (IPTS of them) on the 
SC4020 with the limits of XLF.XRT.YBOT, and YTOP. 

Subroutine LABARC (RI,R,DR,N,NL,KDO) 

LABARC labels to the right-hand side of the 3, 6, 9, and 12 o’clock positions. The value 
is in radial units of specified arcs in the polar grid. 

Input: 

RI The radial offset between the space and the label. 

R The radius in the space of the outermost circle in radial units. 

DR The increment value in radial units to be used for the label. 

N The index to be used for labeling such that every Nth DR will be labeled as 

RI-f-N.DR. 

NL The number of characters to be used in the label. 
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KDO Label Control. 

= 0 Do not Label. 

^ 0 Label according to scheme every Nth DR will be labeled as RI + N.DR. 
OUTPUT: None 

Subroutine LABLV (D,IX,IY,NCHAR,NT,NKMAX) 

LABLV converts a floating point quantity into a number in BCD format with the 
decimal point in the proper position, and it displays the BCD number at the raster 
coordinates specified. 

D The floating point quantity to be printed. 

IX, lY The raster coordinates which will position the first character of the label. This 
may be a leading blank. If the quantity to be displayed is negative, then the 
minus sign will be displayed one character space to the left of (IX, lY). 

NCHAR Number of characters to be displayed including leading blanks and the 
decimal point, if any. NCHAR must be less than 7 or it can be 7 if one of the 
characters is the decimal point. 

NT The number of times each character is to be displayed. NT = 2 will give a 
darker label. 

NDMAX Maximum number of characters to be displayed to left of the decimal point. 

Before the BCD label is displayed, it is right adjusted so that there will be 
NDMAX characters in front of the decimal point. These may now include 
leading blanks. 

NOTE: If NCHAR has a negative sign, then the label will be displayed in scientific 

notation. In this case, NCHAR must be less than 7 and NDMAX can take any 
value. 

If NCHAR = -1 then the display will be of form 
Y.xlO ± YY 

If NCHAR = “3 then the display will be of form 
Y.YYxlO±YY 

If NCHAR = -6 then the display will be of form 
Y.YYYYYxlO±YY 

The space required for scientific notation is greater than for the ordinary rotation. 

Subroutine LABLXC (T,NX,IX) 

LABLXC labels the abscissa axis of a SC4020 plot. 

Subroutine LABLYC (T,NY,1Y) 

LABLYC labels the ordinate axis of a SC4020 plot. 
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Subroutine LINEV (IX1,IY1,IX2,IY2) 


LINEV draws a line between raster coordinates (1X1, lYl) and (1X2, IY2) 

Subroutine MATMUL (A,MA,NA,BM,NA,C) 

MATMUL performs the matrix multiplication C + A*B. 

Subroutine MATPRT (A,M,N,IDENT) 

MATPRT prints a matrix (real or complex) that may have row and column dimensions 
that are unequal. M is number of rows in A, N is number of columns. IDENT is 10 
characters of alphanumeric identification. 

Subroutine MNSQVS(N) 

MNSQVS calculates PSD response, prints frequency and PSD response and plots PSD 
response. MNSQVS also has the entry point FPRINT, 

Subroutine MOMEGA (ROOT,NROOT) 

MOMEGA determines the range of frequencies that will be within the unsteady 
frequency response. 

Subroutine MOVEFR (OMEGAl,DBl,PHIl,OMEGA2,CPS2,FMAG2,DB2,PHI2) 

MOVEFR moves the calculated set of unsteady frequency response data points to a 
higher core address in central memory. This action is performed at the end of the 
QRU,6,3 overlay just before entry into the QRU,6,6 overlay. 

Subroutine MSTORE (A,B»N,M) 

MSTORE moves N elements of A into B. Entry point SWAP swaps elements of A and B 
and entry point ROTSM stores B into A by decrementing. 

Subroutine MULLP (C,N,B,BX,ANS,IERR,RES) 

MULLP finds all roots or a single root of a polynomial with complex coefficients. The 
roots are found by Muller’s method of parabolic iteration with deflation. 

C = an array containing coefficients of the polynomial. 

N = degree of polynomial. 

ANS= complex root array. 


Subroutine MUSER (A,M,N,IMAT) 

MUSER is a portion of the user-supplied FORTRAN interface code. SGCtion 6»8 of 
Volume I provides additional information on its use. 
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Subroutine NODEC (XLOW,XUP,YLOW,YUP,NX,NY) 


NODEC determines the number of decimal places for grid labeling numbers. 

Subroutine NORMAL (U,MATORD, ROW, COL, GROW, GCOL,IPIV) 

NORMAL normalizes the matrix polynomial using Gaussian elimination. This 
subroutine has the same effect as multiplying the matrix polynomial by the inverse of 
the leading matrix. 


Function NXV (X) 

NXV is a function subprogram that converts a floating point abscissa value X to a fixed 
point raster value NXV, ready to plot. A previous call to XSCALV must have been 
made before using IX = NXV(X). If X lies outside the limits, then IX = 0, otherwise, 
NXV = A*X + B. 


FUNCTION NYV (Y) 

NYV is similar to NXV except a fixed point raster is calculated for the ordinate 
direction. 


Subroutine ORDER (ROOT,NROOT) 

Subroutine ORDER orders NROOT eigenvalues, ROOT array, in descending complex 
absolute value order. If the ordering occurs during a COMPUTE FLUTTER option, 
eigenvalues are ordered according to flutter velocity. 

Subroutine PEAKS (N) 


PEAKS calculates the minimums and maximums of an input frequency response. 

Subroutine PGSET (U,MO,MD,PHASE,GAIN) 

PGSET calculates a complex value of gain given a gain and phase value. This complex 
gain value is multiplied times each U matrix element that is designated as a gain 
location. Sections . 3 .3of Volume I has additional detail on root locus. 

Subroutine PLOTTM (TIME,GAIN,N) 

This subroutine plots time responses generated by INLAP onto the SC4020. In addition, 
the plot vectors associated with time response are placed on disk unit 92. 

Subroutine PLOTV (IX,IY,NS) 

This routine plots at point (IX, lY) the character represented by NS. The following table 
shows the possible characters that PLOTV can plot. 
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NS 
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NS 
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NS 
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NS 
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Subroutine PLTLI (N,X,Y,JX,JY,IERR) 

This routine joins N points (X,Y ) by straight lines. If points lie off the grid area, a line will be 
drawn to the intersection(s) of the grid boundary and the line if this intersection(s) exists. 

N = Number of points to be plotted, 

X,Y - Arrays containing the X and Y floating points 
to be plotted. 

JX = X skip factor. 

JY = Y skip factor 

lERR = Error location which, at the end of the program, will contain the number of 
points outside the plotting area. 
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Note that scale factors must be established before using this program. This can be done, 
for example, by a previous call to GRIDIV. 

Subroutine PMRGIN (N) 

PMRGIN calculates the phase margin of an input frequency response. 

Subroutine POINTD (IX, lY, ± NS) 

POINTD plots a symbol, according to the value of NS, at the point specified by (IX, I Y). 
IX,IY = Raster coordinates at the point to be plotted, fixed point. 

±NS = NS is an integer which selects the plotting symbol. If NS is positive an extra 
plotting dot will also be plotted. 

The same table of symbols is used by both POINTV and POINTD. 

Note: This routine should be used whenever possible in place of POINTV(IX,IY,NX,ANY). 


NS Symbol 
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NS Symbol 

NS Symbol 
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NS Symbol 
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Subroutine POINTV (X,Y, + NS) or (IX, lY, INS, ANY) 


POINTV plots a symbol, according to the value of NS, at the point specified by (X,Y) or 
(IX,IY). If floating point arguments are used and the point lies outside the plotting area, no 
symbol is plotted. 

X,Y = Floating point data values of the point to be plotted. Scale factors must 
previously have been established. This can be done, for example by a previous 
use of GRIDIV. 

IX = Raster coordinates of the point to be plotted, fixed point. 

lY = Note:IX 0, lY 0 

±NS = NS is an integer which selects the plotting symbol. If NS is positive, an extra 
plotting dot will also be plotted. 

ANY = This is a dummy argument. 




Note that the table included in POINTD shows the symbols obtained by various values 
of NS 


Function POLAR (RLAI»R1»AA,Y,X) 

POLAR returns the value in the cartesian space of the position of a given polar point. 
POLAR also has the entry point POLER. 

INPUT: 

RI = Radial offset in radial units. 

A1 = Angular offset in degrees. 

Rl = Input radius in radial units. 

AA = Input angle in degrees. 

OUTPUT: 

Y = Vertical value of the point in the cartesian space. 

X = Horizontal value of the point in the cartesian space. 

Y = (Rl - RI) SINE (AA-Al) 

X - (Rl - RI) COSINE (AA - Al) 

Note: Polar does not do any plotting. 

Subroutine PRECON (A,B,D,IA,IB ,IC,ID,NSM,M) 

PRECON scales the matrix A for hessenberg transformation and QR rooting. 

Subroutine PRINTV (N,BCDTXT,IX,IY) 
or (N,nH....,IX JY) 

PRINTV prints N characters in typewriter modestartingat point (IX, lY). Writing will take 
place in horizontal rows, 128 characters in a row and 64 rows to a frame. Care must be taken 
to ensure that the writing does not run off the frame as it will continue at the top left-hand 
corner of the same frame. 

N - Number of characters to print. 

BCDTXT = An array containing the BCD test to be printed, read in with an A type format. 
nH.... = A Hollerith argument containing the text to be printed. 

IX, lY = The raster coordinates for the center of the first character to be printed, 

remembering that each character occupies 8 horizontal raster positions and 16 
vertical raster positions. 
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Subroutine PRUNE (X,Y,M,SYMB) 


PRLINE is part of the printer plot package in QR and is called by PRPLOT, The purpose is to 
scan the X-Y vector and determine if any points would fall onto the current line that is to be 
printed. 


Subroutine PRPLOT (X,Y,M,XP,YP,MP,XZ,YZ,MZ,XL,XR.YT,YB,ITYPE) 

PRPLOT is called when a printer plot is to be generated for time response, root locus, or 
V-g analysis. 


Subroutine PTPLOT (X,Y,N) 

PTPLOT is called when a time response printer plot is to be generated. 

Logical Function QCKROO (GAIN, PHASE) 

QCKROO calculates and roots a polynomial formed from the ZEROS, POLES and gain 
of a transfer function. 

Subroutine QQQQ7 (V,MATORD,MATDEl,VN,MATSIZ,IROW) 

QQQQ7 reforms the V array into VN based on the number of permutations noted in 
array IROW. This subroutine is called after a DELETE AND REDUCE or TRUNCATE 
AND REDUCE input option. 


Subroutine QRCOST 

QRCOST calculates the CRU cost of the QR execution and writes the collected 
information onto a permanent file. 

Logical Function QRROOT (U) 

QRROOT sets up a complete real or complex matrix polynomial eigenvalue rooting. U 
is the starting location of the matrix for which eigenvalues will be found. 


Subroutine RADIIV (A1,D1,R,DR,LRAD) 

RADIIV draws the radii lines of the polar grid at specified angles, with the option of 
labeling. For labeling, RADIIV assumes the space is counterclockwise. It will label as 
though the space is clockwise or counterclockwise from control information. Set the space 
clockwise after the call if the desired space is clockwise. 
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INPUT: 


A1 ~ value of the angle to be displayed at the 3 o’clock position module 360. 

D1 = Increment to the next angle to be labeled or drawn in the counterclockwise 

direction. 

R = Outer-most point of the radii drawn. Distance from the center in radial units. 

DR = Inner-most position for drawing the radii. 

LRAD = Labeling control. 

= 0 No labels will be shown. 

Labels will appear at the outside of every radii. 

Note: If D1 0 no radii will be drawn. 

Subroutine RDMAT (U,MATORD,A,MAT2) 

RDMAT sets up calls to the subroutine READTP 

Subroutine RDVEC (V,MATORD) 

RDVEC sets up calls to the subroutine READTP to read input vectors. 

Subroutine READEM (U,MATORD,MATDEl, VECTOR) 

READEM reads and processes information associated with READ MATRIX TAPE, 
READ VECTOR TAPE, and SPACE TAPE. Entry READFT reads k-value matrices and 
formulates the matrix polynomials corresponding to each individual k-value. Additional 
information on the COMPUTE FLUTTER options can be found in section 6 . 3 . 7of 
Volume I. 

Subroutine REGRUP (U,MATORD,MATDEl,UU,MATSIZ,ROW,COL) 

REGRUP packs the U matrix into the forward part of UU. The rows and columns 
corresponding to ROW(I) or COL(J) equal ".TRUE.” are eliminated. 

Subroutine RESP (OMEGA, FMAG2,DB2,PHI2,CPS2) 

RESP calculates the Db and phase angle of a polynomial ratio (ZEROS(s)/POLES(s)) for 
the input value of omega (radians). If QR is calculating a PSD response, the effects of 
the gust spectrum are included in the calculations. 

Subroutine RHESS (A,IC,IA,N) 

RHESS transforms an upper block triangular matrix into hessenburg from using direct 
reduction with pivoting. The upper portion of A is replaced by the hessenburg form and 
transformation information replaces the lower portion of A. 
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Subroutine RLGRID 


RLGRID reads data associated with the PLOT ROOT LOCUS data card. Section 6,3.3, 
Volume I explains additional detail for this data card. 


Subroutine ROOTS (U,MATORD,MATDEl,II,JJ,ROOT,NR,RAD) 

ROOTS obtains the roots of matrix element II, JJ. The leading coefficient is multiplied 
by the highest power of S. 

Subroutine RPRINT (ROOT, PHA, GAIN, NROOT,NGNCK) 

RPRINT prints the eigenvalues residing in the ROOT array. 

Subroutine RQREIG (A,ROOT,IC,NSM,MM) 

RQREIG computes the eigenvalues of the hessenberg matrix A using the QR algorithm. 
See reference 3 for additional information. 

Subroutine RSTORE (A) 

RSTORE reads the matrix file that has been written by subroutine CSAVE. The 
subroutine is called when the QR program requires that the A matrix be read from disk 
and placed in core. RSTORE also has the entry point RSTOREZ, 

Subroutine RUW (IFOR,IVR,NRWA10) 

RUW is used to read all interactive data in the QR program. 

Subroutine RWAIO (ICARD,N) 

RWAIO reads and writes N alphanumeric CDC words into array ICARD. Entry point 
RAIO reads N words, entry point WAIO writes N words. 

Subroutine RWE12 (A,NCARD,NVARIB,NDIM) 

RWE12 reads and writes variables in format 6E12.5. 

Subroutine RWI5 (N, II, 12, 13, 14, 15, 16) 

RWI5 reads and writes variables in 15 format. 

Subroutine SCALEl (Y,NP,YPMAX,YPMIN) 

SCALE 1 finds linear scaling values that will produce SC4020 plots with convenient X 
or Y axis increments. 


130 



Subroutine SCALE2 (Y,NP,YPMAX,YPMIN) 


SCALE2 finds linear scaling values that will produce printer plots with convenient X or 
Y axis increments. 


Subroutine SCAMB 

SCAMB is part of the interactive error recovery capability of QR. If an error occurs 
while interactive input data is being decoded from alphanumeric to binary, SCAMB is 
entered, an error flag is set, and control is returned to the subroutine RUW, 

Subroutine SENSOR (U,MATORD,MATDEl,PHI,NX) 

SENSOR is entered with the SENSOR and BODY STATIONS control cards. Section 
6 . 3 . 2of Volume I has additional information on these program options. 

Subroutine SETMIV (MTL,MTR,MTB,MTT) 

SETMIV allows a programmer to alter the standard margin assignments used in 
GRIDIV. The standard values are 24 raster counts at the top, left and bottom of the 
grid and 0 at the right. 

MTL Width of area for left margin, in raster counts 
MTR Width of area for right margin, in raster counts 
MTB Width of area for bottom margin, in raster counts 
MTT Width of area for top margin, in raster counts 

This routine makes it possible to provide margin space for multiple lines of printed 
titles and headings. GRIDIV guarantees that these margin values, assigned in SETMIV 
will not be decreased although they may be increased if necessary. 

Subroutine SETORD (U,MAT0RD,MATDE1,10RD,IDJSIZE) 

SETORD determines the apparent state variable order of matrix polynomial U. In 
addition the columns of U are arranged to form the first order state variable matrix. 

Subroutine SORT (X,Y,N) 

This routine sorts the vector X into descending order. If two X elements are 
interchanged, the corresponding Y elements are also interchanged. There are N 
elements. 

Subroutine STUFF (B,C,D,AI,TR,BR,BS,CS,AIS,TRS,AIF,MATORD,MATSIZ, 

MDIFF,MORIG,ISLIDE,GROW) 

STUFF packs the B, C, D, AI, TR, and BR matrices into the front of the U array. 
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Subroutine SUBJOM (SO,UO,Ul,U2,N) 

SUBJOM is used in the unsteady frequency response option. The SO (N x N) matrix is 
formed by substituting a value of frequency, j, into the matrix polynomial U2*S^ + 
U1*S + UO. S is the complex differential operator. 

Subroutine SVFORM (A,MATORD,MATDEl,ID.ROW,COL,GROW,GCOL,IPIV) 

SVFORM transforms the A matrix into state variable form. 

Subroutine TAPDIM (NARRAY,NSIZE,MATDIM) 

TAPDIM is called immediately after a call to subroutine READTP. TAPDIM checks for 
subscript and array length sizes. 

Subroutine TPRINT (N) 

TPRINT controls the detail of time response printing. See sectionS . 3. 5of Volume I for 
additional detail. 

Subroutine TRANS (A,MATORD,MATDEl,VECTOR.LENGTH) 

TRANS is used in the unsteady frequency response option. It transposes the A matrix if 
all of the gains are along one row and builds the vector of knowns for solution in 
lORESP. 

Subroutine UGMARG (OMEGA, GDB, PHI) 

UGMARG is the same as GMARGN except that it is used with the unsteady frequency 
response. 

Subroutine UNSQVS (OMEGA, CPS, FMAG,DB, PHI, PSD, PERCNT) 

UNSQVS is the same as MNSQVS except that is is used with the unsteady frequency 
response. UNSQVS also has the entry point UPSP. 

Subroutine UPEAKS (OMEGA,G, ANGLE) 

UPEAKS is the same as PEAKS except that it is used with the unsteady frequency 
response. 

Subroutine USERSB (A,B,M,N,INSE) 

USERSB is a portion of the user-supplied FORTRAN interface code. Section 6.8 of 
Volume I provides additional information on its use. 
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Subroutine VECTRC (IX1,±IDX,IY1,±IDY) 


VECTRC can be used in place of LINEV when IDX and IDY are less than 64 rasters 
each. Since it will result in significant savings of CPU time, use of VECTRC is 
recommended when there are several short lines to draw. 

Subroutine VPRNTD (N,BCDTXT,IX,IY) or (N,nH.,.,IX,IY) 

VPRNTD does not use the typewriter mode, but prints N characters starting at point 
(IX, lY). Writing will take place in vertical columns. The arguments are the same as 
described for PRINTV, 

Subroutine VREAD (VEC,MATORD,MATDEl) 

VREAD is called when VECTOR or REPLACE VECTOR data cards are being 
processed by QR. Section 6 ,3 . 2of Volume I provides additional detail on these data 
cards. 


Subroutine VUSER (B,M,N,IVEC) 

VUSER is a portion of the user-supplied FORTRAN interface code. Section 6*8 of 
Volume I provides additional information on its use. 

Subroutine WARN (NUM,ID) 

WARN is called when QR has an indication there may be more elements in an array 
than program defaults allow. 


Subroutine WRITER (A,N,K) 

WRITER writes elements of the A array row by row. N rows are written. K represents 
the power of S associated with the A matrix. 

Subroutine WRITUM (U,MATORD,MATDEl) 

WRITUM sets up calls to the subroutine WRITER. Entry point PNCH punches the U 
matrix in QR MATRIX input format. SectionG. 3 .2of Volume I explains the QR PUNCH 
format. 
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XSCALV will compute scale factors for the X coordinate and will save these for use in other 
subroutines. 

XL = Value of X at left-most limit of the scaled plotting area, floating point. 

XR = Value of X at right-most limit of the scaled plotting area, floating point. 

ML = Margin space, in raster counts, to be saved at left of grid. 

MR = Margin space, in raster counts, to be saved at right of grid. 

Only linear scaling is available. 

Subroutine YSCALV (YB,YT,MB,MT) 

YSCALV will compute scale factors for the Y coordinate and will save these for use in 
other subroutines. 

YB = Value of Y at the bottom limit of the scaled plotting area, floating point. 

YT = Value of Y at the top limit of the scaled plotting area, floating point. 

MB = Margin space, in raster counts, to be saved at bottom of grid. 

MT = Margin space, in raster counts, to be saved at top of grid. 

Logical Function ZFIND {U,MATORD,MATDEl) 

ZFIND modifies elements in the U matrix polynomial array so that the matrix when 
rooted will represent the open loop zeros. Sec, 4.2 of Volume I explains the 
automatic formation of poles and zeros. 
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3.0 EXTENT OF CHECKOUT 


The following matrix displays program options used by the QR test cases. 


OPTION 


CASE 


1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 


AIRFORCE COEFFICIENTS 
BEGIN 

BODY STATIONS 

CANCEL 

COMPUTE 

COMPUTE CONVOLUTE 
COMPUTE FLUTTER 

COMPUTE FLUTTER AND ROOT 
LOCUS 

COMPUTE FREQUENCY RESPONSE 
COMPUTE FLUTTER 2 
COMPUTE LOCUS 
COMPUTE POLES 
COMPUTE PSD 
COMPUTE ROOT LOCUS 

COMPUTE ROOT LOCUS WITH 
PREVIOUS POLES 

COMPUTE ROOT LOCUS WITH 
PREVIOUS ZEROS 


X 




X 


X 



X 


X 


COMPUTE ROOT LOCUS WITH 
PREVIOUS POLES AND ZEROS 


COMPUTE TIME RESPONSE 

COMPUTE TIME RESPONSE WITH 
TAU 

COMPUTE UNSTEADY FREQUENCY 
RESPONSE 


X 


X 


X 


X 


COMPUTE ZEROS 
CONTINUOUS 
CRAMER REPLACE 
DC GAIN 
DELETE 

DELETE AND REDUCE 
DIFFERENTIATE 
EXECUTE PLOTS 
EXECUTE PLOTS WITH POLES 
EXECUTE PLOTS WITH ZEROS 


X 


X 


X 


X 


X 

X 


X 

X 

X 
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OPTION 


CASE 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

EXECUTE PLOTS WITH POLES 
AND ZEROS 












X 





FORCING FUNCTION 




X 













FORM POLYNOMIAL 




X 








X 





FORM RATION OF POLYNOMIALS 




X 










X 

X 


FREQUENCY RESPONSE 


X 


X 







X 



X 

X 

X 

GAINS 

GUST SPECTRUM 


X 


X 





X 

X 

X 

X 

X 

X 

X 

X 

INTEGRATE 




X 











X 


MATRIX 

MATRIX COMPLEX 
NO FORM POLYNOMIAL 
PLOT FLUTTER 

X 

X 


X 


X 

X 

X 

X 

X 


X 


X 

X 


PLOT FREQUENCY RESPONSE 


X 









X 

X 


X 

X 

X 

PLOT PSD 




X 











X 


PLOT ROOT LOCUS 

X 



X 







X 

X 


X 


X 

PLOT TIME RESPONSE 




X 











X 

X 

PRINT 

PRINT PARTIAL MATRIX 




X 

X 

X 

X 





X 

X 




PTITLE 




X 








X 




X 

PUNCH 

READ MATRIX TAPE 
READ VECTOR TAPE 

X 















X 

REPLACE 

X 

X 

X 




X 





X 



X 

X 

REPLACE COMPLEX 







X 






X 




REPLACE VECTOR 
SENSOR 













X 



X 

SPACE TAPE 
















X 

STATE 

STOP 

TRUNCATE AND REDUCE 




X 












X 

TIME INTERVALS 




X 











X 

X 

TIME RESPONSE 
TRUNCATE 




X 


X 









X 


TITLE 




X 










X 


X 

UNSTEADY AIRFORCE 
COEFFICIENTS 


X 

X 














VECTOR 

CALLMUSER 

CALLUSERSB 

CALLVUSER 


X 

X 






X 


X 



X 

i 

X 
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